Cold, dilute, trapped bosons as an open quantum system 
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We present a master equation governing the reduced density operator for a single trapped 
mode of a cold, dilute, weakly interacting Bose gas; and we obtain an operator fluctuation- 
dissipation relation in which the Ginzburg-Landau effective potential plays a physically trans- 
parent role. We also identify a decoherence effect that tends to preserve symmetry, even when 
the effective potential has a "Mexican hat" form. 
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The recent observations of Bose condensation in mag- 
netically trapped gases (2|3|,||,[|,|| have brought phenom- 
ena hitherto reserved to condensed matter physics into 
the interdisciplinary field now growing between quantum 
optics and mesoscopic physics. In particular, this break- 
through may offer a weakly interacting system in which 
to study one of the most interesting of quantum phenom- 
ena: the spontaneous breaking of number eigenstates, 
which are invariant (up to a Hilbert space phase) under 
the rotation U = e l8n , into coherent states, which trans- 
form as U\a) — \e l6 a). Such symmetry breaking has 
long been thought to be a basic cause of superfluidity 
and superconductivity. Yet the breaking of the symme- 
try whose Noether charge is particle number is still rather 
more mysterious than the symmetry breaking involved 
in, say, the Standard Model of particle theory, because 
instead of being due to an instability in the microscopic 
Hamiltonian of the system, it is thought to be driven by 
an effective potential of the form 



V= (E 1 -^\ 2 +E 2 \^\ 4 



(1) 



where if) is a second-quantized destruction operator. 

While it has long been customary to include an ef- 
fective potential of this form in the Hamiltonian for a 
gas of weakly interacting bosons {jj , the presence of the 
chemical potential fi must remind us that this effec- 
tive potential (EP) is fundamentally a property of the 
system's thermal environment. In non-equilibrium dy- 
namics, therefore, the EP should most properly appear, 
like temperature, in the fluctuation-dissipation relation 
(FDR). To realize such a result, we must consider Bose 
condensation, and the possibly concomitant symmetry 
breaking, to occur in an open quantum system. The 
open system paradigm is not only natural for gases that 
are cooled by evaporation and thereafter steadily leak out 
of their traps; it is in general a powerful lingua franca 
for describing the rapidly growing common ground be- 
tween quantum optics, mesoscopic quantum mechanics, 
and condensed matter physics. Recent work using this 
paradigm has developed analogs of laser theory for atoms 
P,p|pO|] and excitons [|ll| . The nonlinearity of the inter- 
action between hard spheres makes open system calcula- 
tions difficult for gases in general, but in this Letter we 



present an idealized model of a cold, dilute, trapped Bose 
gas, in which the problem may be solved analytically to 
leading order in small parameters. We are thereby able 
to obtain the fully nonlinear, fully quantum mechanical 
FDR, and so make contact between the open system ap- 
proach and traditional many-body formulations based on 
an EP. 

Instead of the quadratic confining potentials of real 
experiments, we imagine a Bose gas subjected to a deep 
but narrow spherical square potential well, tuned so as 
to possess exactly one single-particle bound state. In the 
second-quantized formalism, this single-particle state be- 
comes a normal mode, with a discrete set of energy levels. 
We treat this mode as an open quantum system, inter- 
acting, via a short-ranged two-particle potential, with a 
reservoir consisting of the continuum of unbound modes. 
We then outline a derivation (to be presented in detail 
elsewhere) that uses Feynman's ordered operator calcu- 
lus |T^] to obtain a Markovian master equation for the 
reduced density operator of the bound mode. This equa- 
tion contains condensate growth and depletion terms 
whose relative strength is characterized by an operator 
fluctuation-dissipation relation, which is physically trans- 
parent, and in which the Ginzburg-Landau EP may be 
discerned. We also find an additional phase diffusion 
term, representing quantum decoherence due to scatter- 
ing of unbound particles by the condensate. We will dis- 
cuss the implications of this term for symmetry breaking 
at the end of this Letter. 

We therefore begin with the second-quantized Hamil- 
tonian 



H = —I dV\Vip f 
2m 



■ - 2\ 2 0(ir - \r)fti> 



(2) 



where 9 is the step function, and a is the scattering 
length of the gas particles (a = lQira 2 being the cross 
section for two hard-sphere bosons scattering with mo- 
mentum transfer small compared to a -1 ). We then di- 
agonalize the quadratic part of H, by defining tp(r) = 
J2kim ^kiraU k i m {r), such that 

(V 2 + k 2 )u k i m = 2\ 2 9(ir- \r)u klm . (3) 
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The Ukim with real k are unbound modes, scattered by 
the well We will let the spectrum of k approach a 
continuum by implicitly assuming a boundary condition 
at large radius. There is also exactly one bound solution, 
which is spherically symmetric; by solving a transcen- 
dental equation we can compute its binding energy to be 
Eb = 0.457 We denote this bound mode wave func- 

tion by ub(t), and its creation and destruction operators 
by and ipB- The bound mode will then constitute the 
observed system, and the continuum modes the environ- 
mental reservoir, which will be in self-equilibrium with 
temperature (fcs/3) -1 and chemical potential p, though 
it may be far from equilibrium with the bound mode. 

We then re-write the Hamiltonian (Q) in terms of the 
normal modes, splitting it into bound mode, continuum, 
and interaction parts: 

H = Hb + He + Hj 
H B = -E b n B + E r h B {n B - 1) 

h 2 

^ C = 2m ^ k2 ^kim^ki™ + H scat = H kin + H scat , (4) 

klm 

where fis = wb^b, and the condensate self- repulsion en- 
ergy E r is of order a?i 2 A 3 TO _1 , proportional to J d 3 ru B . 
H scat is a quartic operator which produces two-particle 
scattering among the continuum modes. It will be con- 
venient to use the interaction picture, in which operators 
evolve under Hb and He- 



[v 2 (k,k',i) ^IJU^Jb^ 



■hk 2 - 



-iAt 



1M0) 



(5) 



Note that A = A(n) is simply the frequency difference 
between the nth and the (re — l)th eigenstate of Hb- The 
operators ^kim(t) evolve under H sca t- This evolution 
is complicated, but for our purposes it turns out to be 
sufficient to note that it makes the two-point function 
decay with time difference. In the case where the gas 
is dilute, we can take sufficient account of this effect by 
approximating 



Skk'Sll'Smm'e 



-7*|t-*'l 



(0) 



where p(/3, p) is the grand canonical ensemble density 
operator, and 7^ = hkcrd/m is the Boltzmann scattering 
rate, for d = ^ (m / 2kTi 2 /3) 3/2 the density of the gas @. 
The equilibration time for the unbound gas is set by the 
thermal average of 7/c, which we will denote by 7. 

In terms of the time-dependent operators, then, we 
can write the interaction picture Hj [using the notation 
k = (fc, I, to)] 

Hiit) = 2^ [&Mlf B kJk 3 + k.c] 



kk'lr, 



(7) 



where all the V n co-efficients are given by integrals of four 
mode functions (uki m an d ub). 

Our goal is to compute the reduced density operator 
of the bound mode, from the initial density operator of 
the bound mode and the reservoir. In the interaction 
picture, this is given by 



p B (t) = Trc 



xlp B (0)®Pc(0W^^ dS " l{s) 



(8) 



where T (T) denote (reverse) time ordering. We have 
assumed that the initial density matrix factorizes, and 
we will take pc(0) to be a grand canonical ensemble 
p(P,/i). In much the same spirit as that of the influ- 
ence functional formalism, we will determine the evolu- 
tion of the reservoir degrees of freedom and perform the 
trace over them before considering the bound mode it- 
self. The bound mode operators ipB(t) are not operators 
in the bath Hilbert space; but we are apparently pre- 
vented from treating them merely as c-numbers in the 
bath Hamiltonian, because they are still operators in the 
bound mode Hilbert space, and their ordering is signif- 
icant. Yet as Feynman observed, this difficulty may be 
overcome by a simple notational trick |T^| . 

In a slight generalization of Feynman's original or- 
dered operator calculus, we here add the device of placing 
primes on all tpB(t) operators that appear to the right of 
the initial density operators in Eqn. (g). We will then 
evaluate the RHS of (g) just as if the bound mode op- 
erators were not operators at all, but keep track of their 
time arguments. We will afterwards be able to restore 
their correct ordering as operators, simply by placing all 
unprimed operators time ordered to the left of ps (0), and 
all primed operators in reverse time order to the right. 
This sort of procedure can always be used, of course; 
but in general it provides only an opaque formal solution 
in terms of time ordered operators. The convenient fea- 
tures of the present idealized model happen to make it 
genuinely powerful here. 

Performing the standard equilibrium analysis after mo- 
mentarily setting a — > confirms that if j3Eb >> 1, Bose 
condensation occurs at very low fugacity ~ e~^ Eb . 
Assuming both /3 and p to be in this regime, then, we can 
evaluate (||) in the dilute gas approximation, where we 
consider only independent two particle collisions. (This 
approximation affects only the unbound bath modes: we 
are not hereby making any assumption concerning the 
bound mode particles, or interactions between them!) 
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We therefore replace Eqn. (m with the dilute gas ex- 
pression 



p B {t) = /5 B (0)exp -- J dsTr[p c (0)(Hi - Hi) 



2/r 



, l , dsTr[p c (0)(Hi - H'j) 



(9) 



da ds'Tr[p c (0)(Hi - ffj) s (ffj - 



o «/o 

where the RHS is to be evaluated using the ordering con- 
vention described above for the bound mode operators. 

The second simplification purchased by our idealiza- 
tions is a time scale separation, which allows us to col- 
lapse the double time integrals in (|J) into single integrals, 
and so obtain Markovian evolution for the bound mode. 
The Boltzman scattering rate is proportional to e^, but 
the condensate grows through collisions between two gas 
particles in the trapping well, which (since the well has a 
finite volume) occur at a rate proportional to e 2/3AI . We 
can therefore use approximations based on the fact that 
the condensate evolves much more slowly than the reser- 
voir equilibrates. Since there are several different terms 
in Hi, we must in fact employ several variations of this 
approximation, and also use the assumption of weak in- 
teraction (Aa << 1). This analysis will be described in 
detail elsewhere, but we will illustrate the least trivial 
and most important step here. 

The following important double time integral from (]p) 
represents a depletion event in which an unbound particle 
of momentum k% dislodges a bound particle, resulting in 
two unbound particles of momenta k 2 and k^ : 



'da fds' e -S(s-s')(fc?-fc!-fc§) e -723| S - S 'l 
ii Jo 



ds + -ip B (s + - e)ipi s (s + - e) 
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= 0(7)+tt / dsip B (s-e)ip j g(s-e) 



2 | ,,2 ,2s A (s + ) + A'(a + )l 



ft H2 ,2 ,2s , A(s) + A'(s) 



^{ — {ki-k 
v 2m v 1 



), (10) 



where 723 denotes jk 2 + 7fc 3 > and in the last line we have 
assumed that t >> 7 _1 . (Nothing significant happens to 
the bound mode prior to this, anyway.) The critical step 
is changing the time arguments of the operators from s 
and s' to s + = (s + s')/2. This can be done because 
when the exponentials in (j^) are expanded, the typical 
time between terms like ( p^| ) is proportional to e~ 2 ^^, but 
the difference in each such term between s and s' can be 
only of order e~^^. Our change of time arguments thus 



only affects the ordering of ips and 4> B operators at order 
e 3/3/i . (Ordering errors involving fis operators must also 
be considered, but this can be shown only to renormalize 

723-) 

The delta function that results from the time scale sep- 
aration enforces self energy {Hb + Hkin) conservation to 
leading order in a 2 , and thus realizes the expectation that 
the dynamics of a dilute gas is dominated by elastic col- 
lisions. (In terms from (^|) with all operators primed or 
all unprimed, there is an imaginary part which is not a 
delta function; but these parts constitute order a 2 A 2 e /3/i 
or smaller renormalizations of the bound mode Hamilto- 
nian, negligible in comparison with the bare Hb-) It is 
crucially important that the delta function we obtain has 
an operator- valued argument. The legitimacy of this cu- 
rious expression, despite its somewhat unnerving appear- 
ance, shows the power of Feynman's ordered operator 
calculus, for it implies the correct operator fluctuation- 
dissipation relation between condensate growth and de- 
pletion. This result seems to be much more difficult to 
obtain by other methods. 

(It is at this point possible to substantially improve our 
dilute gas approximation by taking into account deple- 
tion of the gas modes during condensation. For brevity, 
however, we will here assume that the number of conden- 
sate particles remains a negligible fraction of the number 
of bath particles, so that gas depletion is insignificant.) 

Having in similar fashion collapsed all the double in- 
tegrals in (|^), we obtain for the bound mode reduced 
density operator a Markovian master equation of Bloch- 
Lindblad form. Dropping most of the b subscripts, 
returning to the Schrodinger picture where p is the 
only time-dependent operator, and abandoning the time- 
ordered convention in favour of the ordinary "as- written" 
ordering, it becomes 

■£ = j\P^s\ - a 2 l({fi 2 ,p} - 2npnj 

EtWlQi^}- 2Q±pQl) , (11) 

where a n are constants of order unity which can be eval- 
uated easily to leading order in (ftEb)^ 1 , and the con- 
densate feeding and depleting operators Q± are denned 
to be 



= ,7,1 
■ 



(12) 



This master equation is our main result; it is correct to 
second order in aX and leading order in e^ M . While we 
have made free use of special features of our model to 
derive it, we conclude by making three observations on 
the very general physics it represents. 

Firstly, we note the terms that do not appear in the 
master equation. Hi contains vertices (proportional to 
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Vi and V3) which annihilate one or two bound mode par- 
ticles and create only one or two gas particles (or vice 
versa); but these vertices do not conserve Hb + H^n, 
and so (apart from some more of the negligible correc- 
tions to Hb mentioned above) can only contribute terms 
of order e'^aA) 4 . Excitation of virtual particles from 
the condensate may certainly be less negligible in other 
models; but an energy gap will always tend to suppress 
such processes in favour of elastic ones. 



Secondly it is easy to show that under (11) p relaxes to- 



wards the grand canonical ensemble Z~ 1 e /} ^ n ~ HB \ with 
the full, non-linear Hb- We can consider Eqn. ( |l^ ) to be 
an operator FDR which ensures such relaxation. More- 
over, it is easy to understand this FDR. As we have noted, 
the condensate growth process requires two gas particles 
to collide over the well, and its strength is thus propor- 
tional to e 2/3A \ The depletion process requires only one 
incident gas particle, and so is of order e^, but since 
energy must be conserved to leading order in a 2 A 2 , the 
incoming particle must have high energy in order to dis- 
lodge a bound particle. Depletion is thus enhanced rel- 
ative to growth by a (fugacity) -1 factor, but suppressed 
by a thermal factor. (Since the processes involve two 
Q± operators, the factors in ( p^ ) are the square roots.) 
And our operator delta function in ([To]) implies that 
the thermal factor is precisely e^ A , correctly reflecting 
the fact that interparticle repulsion lowers the binding 
energy as n increases. Furthermore, we can note that 
hA- fj, = V GL (n) - V GL (n - 1), where 

V GL (h) =H-fin= {-E b - + E r ft 2 $ 2 (13) 

is the quantized Ginzburg-Landau effective potential. So 
the Ginzburg-Landau EP simply describes a competition 
between capture of slow gas particles and ejection by fast 
particles, both processes proceeding through elastic col- 
lisions. 

Thirdly, Eqn. ( ^l|) also includes the a-i term, which 
is due to gas particles scattering off the bound mode 
condensate via the elastic part of the V2 vertex in H%. 
This term represents diffusion in the phase canonically 
conjugate to particle number, as may be seen by ex- 
pressing ([ll]) in the Wigner representation, and writing 
q + ip ip = ^/ne ie . The ct2 term then becomes a sec- 
ond derivative with respect to 9. While thinking of the 
real and imaginary parts of ip as independent classical de- 
grees of freedom may suggest that diffusion in 9 should 
be of the form in fact the term does not decay 

with n: It does not represent classical noise, but decoher- 
ence. Yet the argument that a "Mexican hat" form for V 
immediately implies spontaneous breaking of rotational 
symmetry in 9 is actually based on the assumption that, 
when n is large, diffusion in 6 becomes negligible. 

It is certainly possible that decoherence will permit 
spontaneous symmetry breaking in other models (though 
perhaps only for the relative phase(s) between two or 



more condensate modes [[15]]); but scattering of uncon- 
densed particles is a generic process, which may typically 
be expected to produce phase diffusion terms resembling 
that obtained here. These terms may well fail in some 
way to prevent symmetry breaking in other models, but 
they will not obviously vanish due to mere kinematics as 
the number of condensed particles becomes large. This 
strongly suggests the following general conclusion. 

Spontaneous symmetry breaking, in which states of 
definite particle number are replaced by states where the 
phase which is conjugate to particle number is definite 
instead, is important in many areas of physics. A "Mex- 
ican hat" form of effective potential may be a necessary 
condition for such SSB, but it is not sufficient, because 
diffusion associated with decoherence occurs in addition 
to the diffusion characterized by the effective potential. 
This implies that emergence of the symmetry breaking 
order parameter in these cases is not a quasi-classical 
process, but is an especially non-trivial instance of clas- 
sicality emerging from quantum mechanics. 
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